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LITHOGRAPHY SIMULATION METHOD, MASK PATTERN CORRECTION 
METHOD, AND SUBSTRATE TOPOGRAPHY CORRECTION METHOD 

CROSS REFERENCE TO RELATED APPLICATIONS 
5 This application is based upon and claims the benefit of 
priority from the prior Japanese Patent Application No. 
P2003-044443 filed on February 21, 2003, the entire contents 
of which are incorporated herein by reference. 

10 BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention relates to a lithography simulation 
method, a mask pattern correction method, and a processed 
substrate topography correction method, more particularly, 
15 to a lithographic simulation technology for finding the 
distribution of light intensity on the surface of a substrate. 

2. Description of Related Art 

During lithographic processing of a semiconductor device, 
a portion of the light shone onto the surface of the substrate 

20 is sometimes blocked due to uneven portion such as step-like 
parts or other irregularities on the surface of the substrate. 
Conventional lithography simulations calculate the 
distribution of light intensity at the surface of the 
substrate using simple methods that assume the incident light 

25 and exiting light are normal to the surface of the substrate 
(Japanese Laid-Open Patent Publication No. 2002-6475), or 
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using rigorous methods in which the light intensity 
distribution at the step-like part is calculated by solving 
Maxwell's equations. 

The rigorous method is not preferred from a practical 
5 standpoint because the calculations require an enormous 
amount of time. Meanwhile, although the aforementioned 
simple method is able to achieve some degree of calculation 
precision when the shape of the step-like part is somewhat 
gentle , it has poor precision when there are more severe 

10 irregularities on the substrate. The poor calculation 
precision is thought to be caused by the fact that the 
orientation of light reflecting from step-like parts deviates 
from the normal direction, and it has been believed that the 
only way to improve the calculation precision is to use the 

15 aforementioned rigorous method. 

SUMMARY OF THE INVENTION 
An aspect of the present invention provides simulation 
that includes dividing a surface of a substrate onto which 
20 light that is . focused at an aperture angle by a projection 
lens is shone into a first region onto which all of the light 
strikes and a second region onto which a portion of the light 
strikes, calculating an intensity of the light shone onto the 
first region, and calculating an intensity of the light shone 
25 onto the second region . 

Another aspect of the present invention provides a 
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method of mask pattern correction that includes dividing a 
surface of a substrate onto which light that is focused at 
an aperture angle by a projection lens is shone into a first 
region onto which all of the light strikes and a second region 
5 onto which a portion of the light strikes, calculating an 
intensity of the light shone onto the first region, 
calculating an intensity of the light shone onto the second 
region, calculating an exposure pattern to be transferred to 
the substrate using a mask pattern based on the intensities 
10 of the light shone onto the first and second regions, and, 
correcting the shape of mask pattern in such a manner as to 
reduce the amount of deviation of the exposure pattern from 
a desired design pattern. 

Furthermore, another aspect of the present invention 
15 provides a method of mask pattern correction that includes 
dividing a surface of a substrate onto which light that is 
focused at an aperture angle by a projection lens is shone 
into a first region onto which all of the light strikes and 
a second region onto which a portion of the light strikes, 
20 calculating an intensity of the light shone onto the first 
region, calculating an intensity of the light shone onto the 
second region, calculating an exposure pattern to be 
transferred to the substrate based on the intensities of the 
light shone onto the first and second regions, and, correcting 
25 the topology of the substrate in such a manner as to reduce 
the amount of deviation of the exposure pattern from a desired 
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design pattern. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Figure 1 is a block diagram showing a lithography 
5 simulation apparatus for a lithography simulation method 
according to the first embodiment. 

Fig. 2 is a diagrammatic view illustrating the shape of 
a step-like part of the surface of the substrate and the 
projection optical system in a case where the mask pattern 
10 22 is not perpendicular to the planar shape of a polysilicon 
film 28 serving as a base pattern. 

Fig. 3 is an enlarged cross sectional view of the 
step-like part of the substrate 2 9 shown in Fig. 2. 

Fig. 4 is a plan view showing the shape of the pupil 25 
15 and the eclipsed portion 23. 

Fig. 5 is a flowchart of a lithography simulation method 
that uses the lithography simulation apparatus shown in the 
figures . 

Fig. 6 is a plan view showing a layout in which resist 
20 patterns 62 , 63 are arranged such that they are parallel to 
the planar shape of a polysilicon film 61. 

Fig. 7A is a cross sectional view of the second resist 
pattern 63 taken along the line D-D' of Fig. 6* 

Fig. 7B is a cross sectional view of the first resist 
25 pattern 62 taken along the line C-C of Fig. 6. 

Fig. 8 is a diagrammatic view illustrating the shape of 
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a step-like part of the surface of the substrate and the 
projection optical system in a case where the mask pattern 
21 is perpendicular to the planar shape of a polysilicon film 
2B serving as a base pattern. 
5 Fig. 9 is a plan view showing the results of a simulation 

performed using a process simulation method in accordance 
with the second embodiment. 

Fig. 10 is a graph showing the distribution of the light 
intensity in a cross section of the first region 30 indicated 
10 by line B-B' in Figure 9 and the distribution of the light 
intensity in a cross section of the second region 31 indicated 
by line A-A' „ 

Fig. 11 is a cross sectional view illustrating a case 
in which the sloped face 32 of the step-like part is not normal" 
15 to the surface of the semiconductor substrate 27, 

Fig. 12 is a plan view illustrating the refraction of 
the light in a case where light-transmissive materials 41a, 
41b are arranged on the semiconductor substrate 27. 

Fig. 13 is a cross sectional view illustrating the 
20 refraction of the light in a case where a light-transmissive 
material 42 are arranged on the semiconductor substrate 27. 

Fig. 14 is a cross sectional view illustrating a method 
of identifying the eclipsed portion of the light when the 
slope angle of the sloped face of the step-like part varies. 
25 Fi 9- 15 is 3 block diagram of a mask pattern correcting 

apparatus used to implement a mask pattern correction method 
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in accordance with a third embodiment. 

Fig. 16 is a flowchart of a mask pattern correction 
method that uses the mask pattern correcting apparatus shown 
in Fig. 15. 

5 Fig, 17 is a block diagram of a substrate topography 

correcting apparatus used to implement a substrate topography 
correction method in accordance with a fourth embodiment. 

Fig. 18 is a flowchart of a substrate topography 
correction method that uses the substrate topography 
10 correcting, apparatus shown in Fig. 17, 

DETAILED DESCRIPTION OF EMBODIMENTS 
Embodiments of the present invention will now be described 
with reference to the drawings. In descriptions of the 

15 drawings, identical or similar parts will be given identical 
or similar reference symbols. The drawings are merely 
representative examples and do not limit the invention, 
particularly regarding such characteristics as the 
relationship between the thickness and width of individual 

20 layers and the ratios of the thicknesses of one layer to 
another. Furthermore, the invention is not limited by the 
relationships between and/or the ratios of the dimensions of 
one drawing and the dimensions of another. 

9 

(First Embodiment) , 
25 A lithography simulation method in accordance with the first 
embodiment can be achieved using, for example, the 
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lithography simulation apparatus shown in Figure 1. This 
lithography simulation apparatus includes a computing unit 

I having functions for calculating the in-plane intensity 
distribution of light shone onto the surface of the substrate 

5 and at least the following units that are connected to the 
computing unit 1: a projection optical system constant 
storage unit 2, a substrate topography storage unit 3, a mask 
pattern storage unit 4, and a program storage unit 5. The 
computing unit 1 includes a region dividing unit 9 configured 

10 to divide the surface of a substrate - onto which light focused 
to a prescribed aperture angle by a projection lens is shone 
- into a first region onto which all of the light is shone 
and a second region onto which only a portion of the light 
is shone, a first light intensity calculating unit 10 

15 configured to find the intensity of the light shone onto the 
first region, and a second light intensity calculating unit 

II configured to find the intensity of the light shone onto 
the second region. 

The projection optical system constant storage unit 2 stores 
20 projection optical system constants required for lithography 
simulations. The projection optical system constants include 
the pupil function P<5,n) of the optical system, the aperture 
function p 0 ($, r\) of the optical system, the numerical 
aperture NA p of the projection lens, the aperture angle a of 
25 the projection lens, the wavefrpnt aberration w($, t\) of the 
projection lens, the partial coherence ratio (coherency) o 
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(which indicates the ratio between the numerical aperture of 
the illumination optical system and the numerical aperture 
of the projection optical system) , the Fourier spectrum a, 
and the transmission cross-coefficient (TCC) T( fl, f2; gl, 
5 g2) . The substrate topography storage unit 3 stores data 
related to the topography of the substrate. The substrate 
topography includes three-dimensional shape data for the 
substrate, such as the shape of step-like parts resulting from 
irregularities in the surface of the substrate. The mask 

10 pattern storage unit 4 stores data related the shape of the 
mask pattern that will be transferred to the surface of the 
substrate by shining light. The program storage unit 5 stores 
a computer-readable program for calculating the distribution 
of the light intensity at the surface of the substrate. The 

15 program can be stored in either a read-only format or in a 
readable/writeable format. 

The computing unit 1 forms a portion of a central computing 
unit (CPU) in a regular computer system. Although omitted 
from the drawings, the computing unit 1 is connected to a main 

20 storage unit configured to temporarily store programs that 
can be read by a computer for calculating the distribution 
of the light intensity at the surface of the substrate and 
data that will be processed by the computing unit 1 . For the 
projection optical system constant storage unit 2, substrate 

25 topography storage unit 3 , mask pattern storage unit 4, and 
program storage unit 5, it is acceptable to use a 
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semiconductor memory (e.g., a semiconductor ROM or a 
semiconductor RAM) , a magnetic disk device, a magnetic drum 
device, a magnetic tape device, or other auxiliary storage 
device* It is also acceptable to use a portion of the main 
5 storage unit inside the CPU. Through an input/output control 
unit 6, the computing unit 1 is connected to an input device 
7 configured to receive data, commands, and other input from 
an operator and an output device 8 configured to output the 
distribution of the light intensity at the surface of the 

10 substrate. A keyboard, mouse, light pen, or flexible disk 
device can be used as the input device 7. A display device 
or printer can be used as the output device 8. 

Figure 2 is a diagrammatic view illustrating the shape of 
a step-like part of the surface of the substrate and the 

15 projection optical system in a case where the mask pattern 
22 is not perpendicular to the planar shape of a polysilicon 
film 28 serving as a base pattern. As shown in Figure 2, the 
substrate 2 9 includes, for example, a semiconductor substrate 
27 and a polysilicon film 28 deposited onto the surface of 

20 the semiconductor substrate 27. The polysilicon film 28 is 
patterned in a prescribed shape and functions as, for example, 
the gate electrodes of MOS transistors. Consequently, there 
are step-like parts between the outside perimeter of the 
polysilicon film 28 and the semiconductor substrate 27. 

25 Light focused at a prescribed aperture angle a by a 
projection lens is shone onto the surface of the substrate 
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29. When a step-like part exists on the surface of the 
substrate, the light will shine onto the polysilicon film 28 
(which is positioned at the top of the step-like part) and 
a shadow will be cast on the semiconductor substrate 27 (which 
5 is positioned at the bottom of the step-like part) . That is, 
a portion of the light focused at the aperture angle a is 
eclipsed by the polysilicon film 28 at the step-like part- 
In other words, only a portion of the exit pupil 25 of the 
projection lens can be seen from the surface of the 

10 semiconductor substrate 27 in the vicinity of the step-like 
part and the remainder of the exit pupil 25 is eclipsed. This 
eclipsed portion is called the "eclipsed portion" 23. The 
eclipsed portion 23 of the exit pupil 25 does not contribute 
to the image of the mask pattern 22. The amount of light that 

15 is eclipsed depends on the height of the step-like part and 
the distance from a vertical line dropped from the apex of 
the step-like part. Meanwhile, all of the light shines on 
the surface of the polysilicon film 28 and on portions of the 
surface of the semiconductor substrate 27 that are distant 

20 from the step-like part. Thus, the intensity of the light 
striking the surface of the substrate 29 varies depending on 
the whether or not the light is eclipsed and the degree of 
the eclipsing when eclipsing occurs. The surface of the 
substrate 29 includes the surface of the exposed 

25 semiconductor substrate 27 and the surface of the polysilicon 
film 28. The region dividing unit 9 divides the surface of 
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the substrate 29 into a region (first region) onto which all 
of the light strikes and a region (second region) on which 
only a portion of the light strikes and the. remainder is 
eclipsed. 

5 When the mask pattern 22 is not perpendicular to the planar 
shape of the polysilicon film 28 (base pattern) , the eclipsed 
portion 23 of the exit pupil 25 overlaps to different degrees 
with a plurality of diffracted lights generated by the mask 
pattern 22. As shown in Figure 2, when the mask pattern 22 

10 is arranged parallel to the planar shape of the polysilicon 
film 28 (base pattern) , the 0-order light 26b and the +l-order 
light 26a of the mask pattern 22 do not overlap with the 
eclipsed portion 23 but the ^1-order light 2 6c does overlap 
with the eclipsed portion 23. The eclipsed portion 23 of the 

15 exit pupil 25 does not contribute to the image of the mask 
pattern 22. Therefore, the intensity of the light should be 
found based on the shape of the exit pupil 25. 

Figure 3 is an enlarged cross sectional view of the step-like 
part of the substrate 2 9 shown in Figure 2, The region 

20 dividing unit 9 of Figure 1 divides the surface of the 
substrate 29 into a first region 30 and a second region 31 
as follows. As shown in Figure 3, the light is focused onto 
the surface of the substrate 29 with an aperture angle a. The 
height of the step-like part, i.e., the thickness of the 

25 polysilicon film 28, is indicated as "h." Consider a case 
in which the exit pupil 25 is viewed from a point located at 
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a distance s from a vertical line dropped from the apex 45 
of the step-like part. The angle 8 formed between a line 
connecting the apex 45 of the step-like part and the point 
at distance s and a vertical line dropped from the apex 45 
5 of the step-like part is defined by Equation (1) . Since the 
slope angle of the step-like part is vertical, this 
explanation will concern a point located at a distance s from 
the lateral face 32 of the step-like part, 
tan 6 = s/h (1) 

10 When 0 is larger than a, there is no eclipsing of the exit 
pupil 25 at the point located at a distance s from the lateral 
face 32. In other words, the first region 30 is formed by 
points that satisfy the condition 0 > a. Conversely, when 
0 is equal to or smaller than ot, eclipsing of the exit pupil 

15 25 occurs at the point located at a distance s from the lateral 
face 32. In other words, the second region 31 is formed by 
points that satisfy the condition 9 £ a. The top surface of 
the polysilicon film 28 is included in the first region 30 
because there are no step-like parts to block the light. The 

20 region dividing unit 9 determines the first and second regions 
30 , 31 based on the projection optical system constants stored 
in the projection optical system constant storage unit 2 and 
the substrate topography stored in the substrate topography 
storage unit 3. The projection optical system constants 

25 include the numerical aperture NA P of the lens, the aperture 
angle cx, and the partial coherence ratio o. The topography 
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of the substrate includes the height h of the step-like part 
and the slope angle of the lateral face 32. 

The first light intensity calculating unit 10 and second 
light intensity calculating unit 11 of Figure 1 calculate the 
5 intensity of the light shone onto the first region 30 and the 
second region 31. The wavefront aberration is expressed as 
W(£, n) using the coordinates ($, n) of the exit pupil 25 of 
the projection lens. The pupil function P(£, n) of the optical 
system is expressed according to Equation (2) . The pupil 

10 function P(£, ri) is set in consideration of the resolution 
of the light and is 0 in the eclipsed portion 23 of the exit 
pupil 25. In Equation (2), P 0 (£, q) is the aperture function 
and is defined to have a value of 1 inside the exit pupil 25 
and a value of 0 outside the exit pupil 25. The item k is 

15 the wave number and is given by k = 2n/X, where X is the 
wavelength* 

P&r,) = / , 0 (#,7)exp(^(^7 7 )) ( 2 ) 
When an optical system forms the image of a transilluminated 
object, the intensity distribution I (x, y) of the image is 

20 expressed according to Equation (3) . 

OB 

h(*>y)= I ijjT(A,f 2 ;g l ,g 1 Wi>g 1 V<f2.g2h*di(f i -Ak+ifa -g 1 )y'W^Adg x dg t 

to 

(3) 

In Equation (3) , a is the Fourier spectrum of the object and 
a* is the complex conjugate of a. T(fl, f2; gl, g2) expresses 
25 the frequency response of the optical system and is called 
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the transmission cross-coefficient (TCC) . It is expressed 
according to Equation (4) . 

In Equation (4), ? 0 ($, n) is the aperture function and S(£, 
5 r|) is a function expressing the effective light source over 
the pupil. Assuming the brightness distribution of the 
effective light source is uniform, S(^, r\) has the value of 
1 inside a circle of radius a and 0 outside said circle, o 
is the partial coherence ratio. This completes a general 

10 description how imaging of the light is handled. The 
description is disclosed, for example, page 279 of Applied 
Optics 1 (published by Baifukan, 1990) by Tadao Tsuruta. 

The second light intensity calculating unit 11 of Figure 1 
finds the intensity of the light shone onto the second region 

15 31, where the exit pupil 25 is eclipsed. The second region 
31 is shaded when viewed from the exit pupil 25 and the 
intensity of the light shone thereon depends on the height 
h of the step-like part and the distance s from a vertical 
line dropped from the apex 4 5 of the step-like part. 

20 The Figure 4 is a plan view showing the shape of the pupil 
25 and the eclipsed portion 23. When the exit pupil 25 is 
viewed from a point located at a distance s from a vertical 
line dropped from the apex 45 of the step-like part as shown 
in Figure 3, the light is blocked by the step-like part at 

25 the eclipsed portion 23 of the exit pupil 25, which is defined 
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by a circle of radius NA p/ as shown in Figure 4. Therefore, 
when performing image calculations, the intensity of the 
light should be found based on the shape of the exit pupil 
25, which is eclipsed in accordance with the distance s from 
5 the vertical line dropped from the apex 45 of the step-like 
part. The distance d from the center 0 of the exit pupil 25 
to the eclipsed portion 23 and the length 2e of the eclipsed 
portion 23 are expressed according to Equations (5) and (6) 
using the numerical aperture NA P of the projection lens, the 
10 height h of the step-like part, and the distance s from a 
vertical line dropped from the apex 4 5 of the step-like part. 

d = sm0 = s/j[s 2 +h 2 ) (5) 

e = ^NA?-d 2 =^NA?-(s 2 +h 2 ) (6) 

Figure 5 is a flowchart of a lithography simulation method 
15 that uses the lithography simulation apparatus shown in the 
figures. The lithography simulation method will now be 
explained with reference to Figure 5. 

(a) First in step SOI, the region dividing unit 9 shown in 
Figure 1 divides the surface of the substrate 29 - onto which 
20 light that has been focused at a prescribed aperture angle 
ot by the projection lens is shone - into a first region 30 
onto which all of the light strikes and a second region 31 
onto which only a portion of the light strikes. More 
specifically, the region dividing unit 9 divides the surface 
25 into a second region 31 that is shaded when viewed from the 
exit pupil 25 and a first region 30 that is not shaded. The 
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dividing is accomplished using the projection optical system 
constants (including the numerical aperture NA P of the 
projection lens and the coherency o) , the height h of the 
step-like part, and the slope angle of the step-like part. 
5 (b) In step S02, the first light intensity calculating unit 
10 of Figure 1 finds the intensity of the light shone onto 
the first region 30. More specifically, using Equation (3), 
the first light intensity calculating unit 10 finds the 
intensity Io(x, y) of the portion of the light focused at a 

10 prescribed aperture angle a by the projection lens that 
strikes the first region 30. 

(c) Finally in step S03, the second light intensity 
calculating unit 11 of Figure 1 finds the intensity of the 
light shone onto the second region 31. More specifically, 

15 using Equation (3), the second light intensity calculating 
unit 11 finds the intensity I 0 (x, y) of the portion of the light 
focused at a prescribed aperture angle a by the projection 
lens that strikes the second region 31. Thus, the 
distribution of the light intensity at the surface of the 

20 substrate 2 9 can be found using the procedure just described. 
As described previously, the distribution of the light 
intensity at the surface of the substrate 29 is calculated 
based on the shape of the exit pupil 25. Consequently, the 
calculation can be accomplished in less time than is required 

25 by conventional simulation methods using Maxwell' s equations.' 
Additionally, the calculation can be accomplished with higher 
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precision than is achieved by conventional simulation methods 
that assume the incident light and reflected light are 
parallel. In short, the first embodiment can provide a fast 
and highly precise lithography simulation method. Also, when 
5 the light intensity distribution is the distribution of light 
intensity in resist, the topography of the resist can be 
calculated quickly and with high precision. 

The lithography . simulation method just described can be 
expressed as a procedure, i.e., a sequence of processes or 

10 operations. Consequently, this method can be converted into 
a computer program that specifies a plurality of functions 
to be performed by the processor inside a computer system so 
that the method can be executed using a computer system. The 
program can be stored on a recording medium that can be read 

15 by the computer and the aforementioned method can be realized 
in a computerized manner by allowing the computer to read the 
recording medium and execute the program. The recording 
medium can be used as the program storage unit 5 shown in 
Figure 1 or read by the program storage unit 5, and the program 

20 can be used to make the computing unit 1 execute various tasks 
according to a prescribed processing procedure. Feasible 
recording media for storing the program include a memory 
device, a magnetic disk device, an optical disk device, or 
other device that is capable of recording a program. Figure 

25 6 is a plan view showing a layout in which resist patterns 
62, 63 are arranged such that they are parallel to the planar 
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shape of a polysilicon film 61. As shown in Figure 6, the 
a polysilicon film 61 of thickness 200 nm is formed on the 
substrate 60 and a vertical step-like part of height 200 nm 
exists around the outside perimeter of the polysilicon film 
5 61. The first resist pattern 62 is formed adjacent to the 
polysilicon film 61. The second resist pattern 63 is formed 
1 |Ltm away from the outside perimeter of the polysilicon film 
61. The first and second resist patterns 62, 63 are both 500 
nm thick and 250 nm wide. The first resist pattern 62 is 

10 disposed inside the second region 31 and the second resist 
pattern 63 is disposed inside the first region 30. 

Figure 7A is a cross sectional view of the second resist 
pattern 63 taken along the line D-D' of Figure 6. As shown 
in Figure 7A, the shapes of the two lateral faces 64, 65 of 

15 the second resist pattern 63 are horizontally symmetrical. 
Figure 7B is a cross sectional view of the first resist pattern 
62 taken along the line C-C of Figure 6. As shown in Figure 
7B, the shapes of the two lateral faces 66, 67 of the first 
resist pattern 62 are horizontally asymmetrical. The reason 

20 for the asymmetry is that the degree of overlap with the 
eclipsed portion 24 is different for the +l-order light 26a, 
0-order light 26b, and -1-order light 26c shown in Figure 2. 
Also, the width of the first resist pattern 62 is larger at 
the surface of the substrate 60 than the width of the second 

25 resist pattern 63 is because the light intensity is weaker 
in the second region 31 than in the first region 30. 
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(Second Embodiment) 

Similarly to the first embodiment, a lithography simulation 
method in accordance with the second embodiment can be 
achieved using, for example, the lithography simulation 
5 apparatus shown in Figure 1. This lithography simulation 
apparatus includes a computing unit 1 having functions for 
calculating the in-plane intensity distribution of light 
shone onto the surface of the substrate and at least the 
following units that are connected to the computing unit 1: 

10 a projection optical system constant storage unit 2, a 
substrate topography storage unit 3, a mask pattern storage 
unit 4, and a program storage unit 5. The computing unit 1 
includes a region dividing unit 9 configured to divide the 
surface of a substrate - onto which light focused to a 

15 prescribed aperture angle by a projection lens is shone - into 
a first region onto which all of the light is shone and a second 
region onto which only a portion of the light is shone, a first 
light intensity calculating unit 10 configured to find the 
intensity of the light shone onto the first region, and a 

20 second light intensity calculating unit 11 configured to find 
the intensity of the light shone onto the second region. 

Figure 8 is a diagrammatic view illustrating the shape of 
a step-like part of the surface of the substrate and the 
projection optical system in a case where the mask pattern 

25 21 is perpendicular to the planar shape of a polysilicon film 
28 serving as a base pattern. As shown in Figure 8, the 
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substrate 29 includes, for example, a semiconductor substrate 
27 and a polysilicon film 28 deposited onto the surface of 
the semiconductor substrate 27. There are step-like parts 
between the outside perimeter of the polysilicon film 28 and 
5 the semiconductor substrate 27. When the mask pattern 21 is 
perpendicular to the planar shape of the polysilicon film 28 
(base pattern) , the eclipsed portion 23 of the exit pupil 25 
of the projection lens overlaps uniformly with the plurality 
of diffracted lights generated by the mask pattern 21. More 

10 specifically, the +l-order light 26a, 0-order light 26b, and 
-1-order light 26c of the mask pattern 21 overlap to equal 
degrees with the eclipsed portion 23 of the exit pupil 25 of 
the projection lens. When this occurs, it is acceptable to 
calculate the intensity of the 0-order light 26b alone and 

15 it is not necessary to calculate the intensity of the +l-order 
light 26a and -1-order light 26c. In short, the lithography 
simulation can be executed using a simple model in which the 
resolution of the projection lens is not taken into 
consideration. The 0-order light 26b of the mask pattern 21 

20 indicates the shape of the illumination light source itself 
and it is sufficient to consider the manner in which the light 
source is eclipsed. The second embodiment will be described 
for a case in which the slope angle of the step-like part is 
normal to the surface of the semiconductor substrate 27. A 

25 case in which the slope angle of the step-like part is not 
vertical will be discussed later with reference to Figure 11. 
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Similarly to the first embodiment, the region dividing unit 
9 divides rhe surface of the substrate 29 into a first region 
30 where the shape of the light source is not eclipsed and 
a second region 31 where the shape of the light source is 
5 eclipsed. 

The first light intensity calculating unit 10 calculates the 
intensity of the light shone onto the first region 30. The 
light shone onto the first region 30 is the entirety of the 
light focused at a prescribed aperture angle a by the 

10 projection lens. The intensity distribution I (x, y) is 
expressed according to Equation (3) . 

The second light intensity calculating unit 11 finds the 
intensity of the light shone onto the second region 31, where 
the shape, of the light source is eclipsed. The second region 

15 31 is shaded when viewed from the light source and the 
intensity of the light shone thereon depends on the height 
h of the step-like part and the distance s from a vertical 
line dropped from the apex 45 of the step-like part. The 
intensity of the light incident on the second region 31 can 

20 be determined based on the ratio of the surface area of the 
portion of shape of the light source that is eclipsed by the 
step-like part to the surface area of the entire shape of the 
light source. 

When the light source is viewed from a point located at a 
25 distance s from a vertical line dropped from the apex 45 of 
the step-like part as shown in Figure 3, the light is blocked 
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by the step-like part at the eclipsed portion 23 of the light 
source shape, which is defined by a circle of radius NAp, as 
shown in Figure 4 . For the purpose of describing the second 
embodiment, the exit pupil 25 of Figure 4 is assumed to 
5 indicate the shape of the light source. Thus, the radius of 
the circuit is the numerical aperture NAi of the illumination 
lens. The distance d from the center 0 of the light source 
shape to the eclipsed portion 23 and the length 2e of the 
eclipsed portion 23 are expressed according to Equations (5) 
10 and (6) using the numerical aperture NAi of the illumination 
lens, the height h of the step-like part, and the distance 
s from a vertical line dropped from the apex of the step-like 
part. Therefore, the surface area As of the eclipsed portion 
23 is expressed according to Equation (7). 

15 A s =x.NAt.%—e*d (7) 

2.7Z 

The surface area Ai of the light source is expressed 
according to Equation (8) . 
A,-fNA\ (8) 

Therefore, the attenuation r(s, h) of the light quantity 
20 is expressed according to Equation (9) as a function of s and 
h. 

r (s, h) - A 8 /Ai (9) 
Thus, the light intensity Ii(x, y) in the second region 31 
(which is shaded by the step-like part) is expressed according 
25 to Equation (10) using the light intensity I 0 (x, y) of the 
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entire light source given by Equation (3) and the attenuation 
r(s, h) of the light quantity given by Equation (9). 

Ii(-s, y) - I 0 (-s, y) * r ( S , h) (10) 
In this way, the second light intensity calculating unit 11 
5 can find the distribution of the light intensity in the 
vicinity of the step-like part. In other words, the intensity 
of the light shone onto the second region 31 is found based 
. on the ratio of the surface area of the light source shape 
as viewed from the second region 31 to the surface area of 
10 the light source shape as viewed from the first region 30. 
The ratio of the surface area of the light source shape as 
viewed from the second region 31 to the surface area of the 
light source shape as viewed from the first region 30 is 
determined by the height h of the step-like part and the 
15 distance s from a vertical line dropped from the apex of the 
step-like part. Also, when the light intensity distribution 
is the distribution of light intensity in resist, the luminous 
exposure of the resist can be calculated and the topography 
of the resist can be calculated. 

20 

(Experimental Verification) 

Figure 9 is a plan view showing the results of a simulation 
performed using a process simulation method in accordance 
with the second embodiment. Figure 9 shows the in-plane light 
25 intensity for a case involving a substrate made up of a 
semiconductor substrate 27 and a polysilicon film 28 formed 
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on top of the semiconductor substrate in which an exposure 
pattern 4 6 that has uniform line widths and is perpendicular 
to the polysilicon film 28 is transferred to the substrate. 
The thickness of the polysilicon film 28, i.e., the height 
5 of the step-like part, is 200 nm. The exposure pattern has 
a line width of 0.25 and is a positive type exposure pattern, 
i.e., a pattern formed where the light is not shone. The 
wavelength of the light is 248 nm, the numerical aperture NA P 
of the projection lens is 0.6, and the coherency o is 0.75. 

10 There are step-like parts between the outside perimeter of 
the polysilicon film 28 and the semiconductor substrate 27. 
Portions of the surface of the semiconductor substrate 27 in 
the vicinity of the step-like part correspond to the second 
region 31 where the light source is eclipsed; all other 

15 portions of the surface of the semiconductor substrate 27 and 
the surface of the polysilicon film 28 correspond to the first 
region 30 where all of the light shines. As shown in Figure 
9, the light intensity is weaker in the second region 31 than 
in the first region 30. As a result, the line width D2 of 

20 the exposure pattern 4 6 in the second region 31 is wider than 
the line width Dl of the same in the first region 30. 

Figure 10 is a graph showing the distribution of the light 
intensity in a cross section of the first region 30 indicated 
by line B-B in Figure 9 and the distribution of the light 

25 intensity in a cross section of the second region 31 indicated 
by line A-A. As shown in Figure 10, in the cross section B-B 
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of the first region 30 shown in Figure 9, intense light is 
shone everywhere except in the domain of the exposure pattern 
46 while the intensity of the light decreases rapidly in the 
domain of the exposure pattern 46. Meanwhile, in the cross 
5 section A-A' of the second region 31 shown in Figure 9, 
intensity of the light shone everywhere except in the domain 
of the exposure pattern 4 6 is approximately half that of the 
cross section B-B' and the intensity of the light decreases 
gradually in the domain of the exposure pattern 46. As a 

10 result, assuming that regions where the light intensity is 
0.3 or less form the exposure pattern 46, the line width D2 
in the second region 31 is wider than the line width Dl in 
the first region 30, These simulation results confirm 
experimentally, based on an experiment conducted by the 

15 inventors, that a process simulation method in accordance 
with the second embodiment faithfully reproduces an actual 
photolithography process. 

(First Variation) 
20 The first and second embodiments have been described for 
cases in which the sloped face 32 of the step-like part is 
normal to the surface of the semiconductor substrate 27, as 
shown in Figure 3 and 8 . However, the invention is not limited 
to such cases and can also be applied to cases in which the 
25 sloped face 32 is not normal to the substrate surface. 

Figure 11 is a cross sectional view illustrating a case in 
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which the sloped face 32 of the step-like part is not normal 
to the surface of the semiconductor substrate 27. As shown 
in Figure 11, light is focused at an aperture angle a onto 
a step-like part of a substrate 2 9 that is made up of a 
5 semiconductor substrate 27 and a polysilicon film 28 
deposited onto the surface of the semiconductor substrate 27. 
The sloped face 32 of the polysilicon film 28 is sloped at 
a prescribed slope angle j5 with respect to a vertical line 
dropped from the apex 45 of the step-like part. When the slope 

10 angle p is larger than the aperture angle a of the projection 
lens, no portion of the exit pupil is eclipsed by the step-like 
part on the surface of the substrate 29 and a second region 
31 is not formed on the surface of the semiconductor substrate 
27 in the vicinity of the step-like part. Conversely, when 

15 the slope angle (J is equal to or smaller than the aperture 
angle a of the projection lens, a portion of the exit pupil 
is eclipsed by the step-like part on the surface of the 
substrate 29 as shown in Figure 11 and a second region 31 is 
formed which is shaded because not all of the light is shone 

20 thereon. In this variation, too, similarly to the first and 
second embodiments, the region dividing unit 9 (shown in 
Figure 1) can divide the surface into first and second regions 
30, 31 and the second light intensity calculating unit 11 can 
calculate the intensity of the light shone onto the second 

25 region 31 based on the height h of the step-like part and the 
distance s from a vertical line dropped from the apex 45 of 
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the step-like part. Therefore, a quick and highly precise 
lithography simulation can be executed even when the sloped 
face 32 is not vertical. 

5 (Second Variation) 

Lithography simulation methods in accordance with the first 
and second embodiments can also be applied to cases in which 
a material that transmits light is coated or deposited onto 
a semiconductor substrate having a step-like part. 

10 Figure 12 is a plan view illustrating the refraction of the 
light in a case where light-transmissive materials 41a, 41b 
are arranged on the semiconductor substrate 27. As shown in 
Figure 12, the substrate includes a semiconductor substrate 
27, a polysilicon film 28 deposited onto the surface of the 

15 semiconductor substrate 27, and light-transmissive materials 
41a, 41b. The light-transmissive materials 41a, 41b are so 
arranged as to contact the lateral faces of the polysilicon 
film 28. The refractive index of the light-transmissive 
materials 41a, 41b is nl and the top surface and lateral faces 

20 of the light-transmissive materials 41a, 41b are parallel and 
normal, respectively, to the surface of the semiconductor 
substrate 27. A portion of the light shone onto the top 
surface of the light-transmissive material 41a at an incident 
angle a enters the light-transmissive material 41a at a 

25 refraction angle (31. A portion of the light striking the 
lateral surface of the light-transmissive material 41a at an 
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incident angle (32 exits the light-transmissive material 41a 
at a refraction angle (33. Thereafter, the light strikes the 
surface of the semiconductor substrate 27 at an incident angle 
of ©1. The angles a, (31, (32/ (33, and col satisfy Equations 
5 (11) and (12) . Thus, the incident angle col at which the light 
strikes the semiconductor substrate 27 can be found using 
Snell's law while considering the topography of the substrate 
and the refractive indices of the light-transmissive 
materials 41a, 41b, making it possible to implement a 
10 lithography simulation method in accordance with the first 
embodiment • 

sin a / sin (3i = sin p 3 / sin (3 2 = ni (11) 
(3i 2 + (3 2 2 - (3 3 2 + Q! 2 - (n/2) 2 (12) 
Figure 13 is a cross sectional view illustrating the 

15 refraction of the light in a case where a light-transmissive 
material 42 are arranged on the semiconductor substrate 27. 
As shown in Figure 13, the substrate includes a semiconductor 
substrate 27, a polysilicon film 28 deposited on top of the 
semiconductor substrate 27, and a light-transmissive 

20 material 42 arranged on top of the polysilicon film 28, The 
light-transmissive materials 42 is so arranged as to contact 
the top surface of the polysilicon film 28. The refractive 
index of the light-transmissive material 42 is n 2 and the top 
surface and lateral faces of the light-transmissive material 

25 42 are parallel and normal, respectively, to the surface of 
the semiconductor substrate 27. A portion of the light shone 
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onto the surface of the light-transmissive material 42 at an 
incident angle oc enters the light-transmissive material 42 
at a refraction angle £i. A portion of the light striking the 
lateral surface of the light-transmissive material 42 at an 
5 incident angle (5 2 exits the light-transmissive material 42 at 
a refraction angle 0 3 . Thereafter , the light strikes the 
surface of the semiconductor substrate 27 at an incident angle 
of g> 2 . The angles a, Pi# 1*2/ p3, and u 2 satisfy Equations (13) 
and (14) . 

10 sin a / sin (Ji ■ sin (3 3 / sin (5 2 = n 2 (13) 

Pi 1 + P2 2 - P3 2 + co 2 2 = (n/2) 2 (14) 
Thus, the incident angles coi and ca 2 at which the light strikes 
the semiconductor substrate 27 can be found using Snell's law 
while considering the topography of the substrate and the 

15 refractive indices of the light-transmissive materials 41a, 
41b, 42, making it possible to execute a quick and highly 
precise lithography simulation. Examples of 

light-transmissive materials 41a, 41b, 42 include oxide films 
(Si0 2 ) and nitride films (Si 3 N 4 ) . 

20 (Third Variation) 

Figure 14 is a cross sectional view illustrating a method 
of identifying the eclipsed portion of the light when the 
slope angle of the sloped face of the step-like part varies. 
As shown in Figure 11, when the slope angle (S of the sloped 

25 face 32 is uniform, the region dividing unit 9 (shown-in Figure 
1) can determine if a second region 31 exists by comparing 
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the slope angle |J with the aperture angle a of the projection 
lens. When the slope angle of the sloped face 4 9 of the 
step-like part varies within the sloped face 4 9 shown in 
Figure 14, however, the region dividing unit 9 must determine 
5 the second region 31 while considering the slope angle at 
every point along the sloped face 49. 

The substrate includes a semiconductor substrate 27 and a 
deposited material 43 deposited onto the semiconductor 
substrate 27. The exposed surface of the deposited material 

10 43 begins to slope gently with respect to the surface of the 
semiconductor substrate 27 at a slope start point 47 and the 
slope angle increases such that said exposed surface is 
approximately normal to the surface of semiconductor 
substrate 27 when it reaches that surface. When the slope 

15 angle of the sloped surface 49 of the deposited material 43 
is iarger than incident angle a of the light, a second region 
where the step-like part casts a shadow is not formed. 
Conversely, when the slope angle of the sloped surface 49 of 
the deposited material 43 is equal to or smaller than incident 

20 angle oc of the light, a second region where the step-like part 
casts a shadow is formed Therefore, the region dividing unit 
9 (shown in Figure 1) compares the slope angle at each point 
on the sloped face with the incident angle a of the light and 
divides the substrate surface into first and second regions 

25 30, 31 based on the distance s from a vertical line dropped 
from each point on the sloped face. In the case shown in Figure 
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14, the point on the sloped face 4 9 where the slope angle of 
the sloped face 4 9 equals the incident angle a of the light 
is defined as the apex 4 8 and the substrate surface is divided 
into first and second regions 30, 31 based on the height h 
5 of the apex 48 and the distance from a vertical line dropped 
from the apex 48. 

Thus, a quick and highly precise lithography simulation can 
be executed even when the slope angle of the sloped face 4 9 
of the step-like part varies within the sloped face 49. 

10 (Third Embodiment) 

Figure 15 is, a block diagram of a mask pattern correcting 
apparatus used to implement a mask pattern correction method 
in accordance with a third embodiment. This mask pattern 
correcting apparatus includes a computing unit 1 having 

15 functions for calculating the in-plane intensity 
distribution of light shone onto the surface of the substrate 
29 and at least the following units that are connected to the 
computing unit 1: a projection optical system constant 
storage unit 2, a substrate topography storage unit 3, a mask 

20 pattern storage unit 4, and a program storage unit 5. The 
computing unit 1 includes a region dividing unit 9 configured 
to divide the surface of a substrate 29 - onto which light 
focused to a prescribed aperture angle a by a projection lens 
is shone - into a first region 30 onto which all of the light 

25 is shone and a second region 31 onto which only a portion of 
the light is shone, a first light intensity calculating unit 
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10 configured to find the intensity of zhe light shone onto 
the first region 30, a second light intensity calculating unit 

11 configured to find the intensity of the light shone onto 
the second region 31, a projection pattern calculating unit 

5 12 configured to calculate an exposure pattern 4 6 to be 
transferred to the substrate 29 using a mask pattern 21 based 
on the intensities of the light shone onto the first and second 
regions 30, 31, and a mask pattern correcting unit 13 
configured to correct the shape of the mask pattern in such 

10 a manner as to reduce the amount of variation between the 
exposure pattern 4 6 and the desired design pattern. 

Through an input /output control unit 6, the computing unit 
1 is connected to an input device 7 configured to receive data, 
commands, and other input from an operator and an output 

15 device 8 configured to output the distribution of the light 
intensity at the surface of the substrate. The program 
storage unit 5 stores in a fixed manner a computer-readable 
program for correcting the mask pattern. 

Based on the light intensity distribution, the exposure 

20 pattern calculating unit 12 calculates the shape of the 
exposure pattern that will be formed on the resist coated onto- 
the substrate 29. The mask pattern correcting unit 13 
compares the desired design pattern, i.e., the pattern one 
wishes to form on the resist, to the calculated exposure 

25 pattern, extracts the amount of variation between the two, 
and corrects the mask pattern 21 in such a manner as to reduce 
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the amount of variation. 

Figure 16 is a flowchart of a mask pattern correction method 
that uses the mask pattern correcting apparatus shown in 
Figure 15. A mask pattern correction method that uses the 
5 mask pattern correcting apparatus shown in Figure 15 will now 
be described. 

(a) First, in step SOI, the region dividing unit 9 shown in 
Figure 15 divides the surface of the substrate 29 - onto which 
light that has been focused at a prescribed aperture angle 

10 a by the projection lens is shone - into a first region 30 
onto which all of the light strikes and a second region 31 
onto which only a portion of the light strikes. More 
specifically, the region dividing unit 9 divides the surface 
into a second region 31 that is shaded when viewed from the 

15 exit pupil and a first region 30 that is not shaded. The 
dividing is accomplished using the projection optical system 
constants (including the numerical aperture NA P of the 
projection lens and the coherency a) , the height h of the 
step-like part, and the slope angle of the step-like part. 

20 (b) In step S02, the first light intensity calculating unit 
10 of Figure 15 finds the intensity of the light shone onto 
the first region 30. More specifically, using Equation (3), 
the first light intensity calculating unit 10 finds the 
intensity I 0 (x, y) of the entirety of the light focused at a 

25 prescribed aperture angle a by the projection lens. 

(c) In step S03, the second light intensity calculating unit 
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11 of Figure 15 finds the intensity of the light shone onto 
the second region 31. More specif ically, using Equation (3) , 
the second light intensity calculating unit 11 finds the 
intensity Io(x, y) of the portion of the light focused at a 
5 prescribed aperture angle a by the projection lens that 
strikes the second region 31. Thus, the intensity 
distribution of the light at the surface of the substrate 29 
can be found using the steps just described. 

(d) In step S04, the exposure pattern calculating unit 12 
10 of Figure 15 calculates the exposure pattern 4 6 to be 

transferred to the substrate 29 using a mask pattern 21 based 
on the intensity of the light shone onto the first and second 
regions 30, 31. 

(e) Finally, in step S05, the mask pattern correcting unit 
15 13 of Figure 15 corrects the shape of the mask pattern in such 

a manner as to reduce the amount of variation between the 
exposure pattern 4 6 and the desired design pattern. 

Thus, as previously explained, the exposure pattern is 
calculated based on the intensity of the light shone on the 

20 first and second regions 30, 31, the exposure pattern is 
compared to the design pattern, and corrections are applied 
to the mask pattern. Consequently, the mask pattern can be 
corrected in less time than is required with a conventional 
simulation method using Maxwell's equations. Furthermore, 

25 the mask pattern can be corrected with higher precision that 
is achieved with conventional simulation methods that assume 
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parallel light. In short, the third embodiment can provide 
a fast and highly precise mask pattern correction method. 
(Fourth Embodiment) 

Figure 17 is a block diagram of a substrate topography 
5 correcting apparatus used to implement a substrate topography 
correction method in accordance with a fourth embodiment- A 
substrate topography correction method in accordance with the 
fourth embodiment can be achieved using, for example, the 
substrate topography correcting apparatus shown in Figure 17. 

10 The substrate topography correcting apparatus includes a 
computing unit 1 having functions for calculating the 
in-plane intensity distribution of light shone onto the 
surface of the substrate 29 and at least the following units 
that are connected to the computing unit 1: a projection 

15 optical system constant storage unit 2, a substrate 
topography storage unit 3, a mask pattern storage unit 4, and 
a program storage unit 5. The computing unit 1 includes a 
region dividing unit 9 configured to divide the surface of 
a substrate 29 - onto which light focused to a prescribed 

20 aperture angle a by a projection lens 25 [sic] is shone - into 
a first region 30 onto which all of the light is shone and 
a second region 31 onto which only a portion of the light is 
shone, a first light intensity calculating unit 10 configured 
to find the intensity of the light shone onto the first region 

25 30, a second light intensity calculating unit 11 configured 
to find the intensity of the light shone onto the second region 
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31, a projection pattern calculating unit 12 configured to 
calculate an exposure pattern 4 6 to be transferred zo the 
substrate 29 based on the intensities of the light shone onto 
the first and second regions 30, 31, and a substrate 
5 topography correcting unit 14 configured to correct the 
topography of the substrate 29 in such a manner as to reduce 
the amount of variation between the exposure pattern 4 6 and 
the desired design pattern. 

Through an input/output control unit 6, the computing unit 

10 1 is connected to an input device 7 configured to receive data, 
commands, and other input from an operator and an output 
device 8 configured to output the distribution of the light 
intensity at the surface of the substrate. The program 
storage unit 5 stores in a computer-readable program for 

15 correcting the topography of the substrate 29. 

Based on the light intensity distribution, the exposure 
pattern calculating unit 12 calculates the shape of the 
exposure pattern that will be formed on the resist coated onto 
the substrate 29. The substrate topography correcting unit 

20 14 compares the desired design pattern, i.e., the pattern one 
wishes to form on the resist, to the calculated exposure 
pattern, extracts the amount of variation between the two, 
and corrects the topography of the substrate 29 in such a 
manner as to reduce the amount of variation. 

25 Figure 18 is a flowchart of a substrate topography 
correction method that uses the substrate topography 
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correcting apparatus shown in Figure 17, A substrate 
topography correction method that uses the substrate 
topography correcting apparatus shown in Figure 17 will now 
be described with reference to Figure 18, 
5 (a) First, in step S01 f the region dividing unit 9 shown in 
Figure 17 divides the surface of the substrate 29 - onto which 
light that has been focused at a prescribed aperture angle 
a by the projection lens is shone - into a first region 30 
onto which all of the light strikes and a second region 31 

10 onto which only a portion of the light strikes. More 
specifically, the region dividing unit 9 divides the surface 
into a second region 31 that is shaded when viewed from the 
exit pupil and a first region 30 that is not shaded. The 
dividing is accomplished using the projection optical system 

15 constants (including the numerical aperture NA P of the 
projection lens and the coherency o) , the height h of the 
step-like part, and the slope angle of the step-like part, 

(b) In step S02, the first light intensity calculating unit 
10 of Figure 17 finds the intensity of the light shone onto 

20 the first region 30. More specifically, using Equation (3), 
the first light intensity calculating unit 10 finds the 
intensity I 0 (x, y) of the entirety of the light focused at a 
prescribed aperture angle Of by the projection lens. 

(c) In step S03, the second light intensity calculating unit 
25 11 of Figure 17 finds the intensity of the light shone onto 

the second region 31. More specifically, using Equation (3) , 
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the second light intensity calculating unit 11 finds the 
intensity Iq(x, y) of the portion of the light focused at a 
prescribed aperture angle a by the projection lens that 
strikes the second region 31. Thus, the intensity 
5 distribution of the light at the surface of the substrate 29 
can be found using the steps just described, 

(d) In step S04, the exposure pattern calculating unit 12 
of Figure 17 calculates the exposure pattern 46 to be 
transferred to the substrate 29 based on the intensity of the 

10 light shone onto the first and second regions 30, 31. 

(e) Finally, in step SOS, the substrate topography 
correcting unit 14 of Figure 17 corrects the topography of 
the substrate 29 in such a manner as to reduce the amount of 
variation between the exposure pattern 46 and the desired 

15 design pattern. 

Thus, as previously explained, the exposure pattern is 
calculated based on the intensity of the light shone on the 
first and second regions 30, 31, the exposure pattern is 
compared to the design pattern, and corrections are applied 

20 to the substrate topography. Consequently, the substrate 
topography can be corrected in less time than is required with 
a conventional simulation method using Maxwell's equations. 
Furthermore, the substrate topography can be corrected with 
higher precision that is achieved with conventional 

25 simulation methods that assume parallel light. In short, the 
fourth embodiment provides a fast and highly precise 
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substrate topography correction method. 

Although the present invention has been described 
heretofore using the first through fourth embodiments and 
variations thereof, descriptions and drawings forming a 
5 portion of this disclosure should not be construed to limit 
this invention. It is believed that this disclosure will make 
various alternative embodiments, working examples, and 
applied technologies obvious to those skilled in the art. 
As described previously, the present invention makes it 
10 possible to provide, a lithography simulation method, a mask 
pattern correction method, and a substrate topography 
correction method, each of which is fast and highly precise. 
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